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As long as magnetic fields remain frozen into the gas, the magnetic braking prevents the forma¬ 
tion of protosteiiar discs. This condition is subordinate to the ionisation fraction characterising 
the inmost parts of a collapsing cloud. The ionisation level is established by the number and the 
energy of the cosmic rays able to reach these regions. Adopting the method developed in our pre¬ 
vious studies, we computed how cosmic rays are attenuated as a function of column density and 
magnetic field strength. We applied our formalism to low- and high-mass star formation models 
obtained by numerical simulations of gravitational collapse that include rotation and turbulence. 
In general, we found that the decoupling between gas and magnetic fields, condition allowing 
the collapse to go ahead, occurs only when the cosmic-ray attenuation is taken into account with 
respect to a calculation in which the cosmic-ray ionisation rate is kept constant. We also found 
that the extent of the decoupling zone also depends on the dust grain size distribution and is larger 
if large grains (of radius ~ 10“^ cm) are formed by compression and coagulation during cloud 
collapse. The decoupling region disappears for the high-mass case due to magnetic field diffusion 
that is caused by turbulence and that is not included in the low-mass models. We infer that a 
simultaneous study of the cosmic-ray propagation during the cloud’s collapse may lead to values 
of the gas resistivity in the innermost few hundred AU around a forming protostar that is higher 
than generally assumed. 
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1. Introduction 


Cosmic rays (hereafter CRs) have a twofold role in the interstellar medium. In fact, they 
both originate the chemistry in molecular clouds and have a fundamental function in controlling 
the collapse of a cloud and the resulting star formation. In prestellar and Class 0 sources, CRs 
represent the main ionising agent, since X-rays ionisation arises only in presence of embedded 
young stellar objects [|T^, |3^ and interstellar UV photons are absorbed for a visual extinction 
Ay > 4 magnitudes [14|. During the collapse of a cloud, magnetic fields coupled to the gas cause 


a brake of any rotational motions, at least as long as the field remains frozen info fhe gas and fhe 
rofafion axis of fhe cloud is close fo fhe mean direcfion of fhe field, e.g. [IHM- This is why fhe 
sfudy of fhe formafion of circumsfellar discs sfill presenfs fheorefical challenges. However, discs 
around Classes I and II young sfellar objecfs are commonly observed [ 0^ , |^], and fhere is also 
some evidence of discs around Class 0 objecfs [ 


18]. 


In order fo alleviafe fhe magnefic braking, differenf mechanisms have been infroduced: (/) 
non-ideal magnefohydrodynamic (MHD) effecls ® iSSil; (//) misalignmenf befween fhe 
main magnefic field direcfion and fhe rofafion axis [^, |^; (///) furbulenf diffusion of fhe magnefic 
field [^, (/v) flux redisfribufion driven by fhe inferchange insfabilify [jlTI]; and (v) depletion 


of fhe infalling envelope anchoring fhe magnefic field [ jl^ p^ ]. 

Non-ideal MHD effecls, namely ambipolar. Hall, and Ohmic diffusion, depend on fhe abun¬ 
dances of charged species as well as on fheir mass and charge. The ionisafion fraction, in furn, is 
determined by CRs in cloud regions of relafively high column density where sfar formafion lakes 
place. The CR ionisafion rale is usually assumed fo be equal fo a “slandai'd” (conslanl) value 
of s“^ 1^]. However, CRs interacting wilh H 2 in a molecular cloud lose energy by sev¬ 


eral processes, mainly by ionisation losses (see [|20|], hereafler PGG09). As a consequence, while 
low-energy CRs {E < 100 MeV) are possibly prevented from entering a molecular cloud because 
of slreaming insfabilify [||, [Tt]], high-energy CRs are slowed down fo energies lhal are relevanl 
for ionisafion (ionisation cross secfions for profons and elecfrons colliding wifh H 2 peak al abouf 
10 keV and 0.1 keV, respectively). 

If has been shown in PGG09 fhal can decrease by abouf fwo orders of magnifude from dif¬ 
fuse regions of column densities ~ 10^' cm“^ (^^2 10“'^ - 10“^^ s“^) fo dense cores and massive 


profosfellar envelopes wifh column densities of ~ 10"^'* cm“"^ (f- ~ 10~^® - 10~^' s“^). A furfher 
affenuafion is caused by fhe presence of magnetic fields. In facf, a poloidal field fhreading a molec¬ 
ular cloud core reduces on average by a facfor 3-4 depending on fhe magnefisafion degree and 
fhe position inside fhe core [00]- A sfronger reducfion of fakes place in fhe inner region of 
a core, where fhe formation of a profoslellar disc is expected fo occur (< 100 - 200 AU) because 


fhe foroidal field componenf generated by rofafion boosfs fhe pafh lengfh of particles [24]. 

In order fo properly freaf fhe problem of fhe influence of CRs on fhe collapse dynamics and 
fhe circumsfellar disc formation, one should compufe fhe CR propagafion self-consislenlly wilh 
fhe evolution of magnetic field and densify, following al fhe same time fhe formafion and fhe de- 
slruclion of chemical species. Since Ihis approach would be exlremely time-consuming from a 
numerical poinl of view, (/) we lake some snapshols of densify and magnefic field configuralion 
from ideal MHD simulalions lhal do nol include any gas resistivity; (//) we follow fhe propaga- 
lion of CRs, compuling fheir spatial dislribulion and fhe CR ionisation rale; (Hi) using a simplified 
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chemical model, we approximately calculate the chemical composition as a function of the values 
of computed in the previous step; (/v) finally, we compute the microscopic resistivities (ambipo- 
lar. Hall, and Ohmic) and compare the time scale of magnetic field diffusion ts to the dynamical 
time scale tdyn at each point in the model to determine the region of magnetic decoupling, where 
tg < tdyn- The hypothesis of ideal MHD on which the simulations are based becomes invalid by def¬ 
inition inside the decoupling region. However, even if our calculations are not fully self-consistent, 
we show that assuming the proper value of at different depths, gives strong constraints on the 
disc formation. 


2. Energy loss processes and magnetic effects on cosmic-ray propagation 


While passing through a molecular cloud, a CR undergoes collisions with molecular hydrogen. 
According to its energy and composition, it is slowed down due to processes that are specific of a 
parficular kind of particle (bremssfrahlung, synchrofron emission, and inverse Compton scattering 
for elecfrons; elasfic inferacfions, pion producfion, and spallafion for profons) or common bofh 
to CR profons and elecfrons (Coulomb and inelastic inferacfions, and ionisation). PGG09 show 
fhaf even if a local CR infersfellar specfrum is lacking of low-energy particles, fhe slowing-down 
of high-energy CR profons and elecfrons during fheir propagafion produces a low-energy fail. Our 
modelling is able to explain fhe decrease of wifh increasing hydrogen column density compufed 
from observations. In particular, a proton componenf af low energies, and mosf likely also an 
elecfron componenf, could be necessary fo reproduce fhe dafa. 

The quanfify fhaf describes fhe energy losses during fhe propagafion is called energy loss 
function and if is given by 


LkiEi,) - - 


dA(H2) ’ 


(2.1) 


where A^(H 2 ) is fhe column density of fhe medium in which fhe particle of species k and energy E^ 
propagafes (see e.g. leff panel of Fig. 3 in [p^). Besides energy losses, one has to accounf for fhe 
facf fhaf CRs are charged parficles and fhey move along fhe field lines following an helicoidal pafh. 
This means fhaf fhey “see” a larger H 2 column densify wifh respecf fo a rectilinear propagafion, so 
fhaf fhe column densify in Eq. (pl|) reads 


N{a) = 



n{e)di, 


( 2 . 2 ) 


where f’^ax is the maximum depth reached inside the core and n{i) is the H 2 volume density. The 
angle a, called pitch angle, is the angle between the CR velocity and the direction of the magnetic 
field. Its evolution during the CR propagation reads 


a = arccos 


y/i “T+TCOs^aicM, 


(2.3) 


where x is the ratio between the local and the interstellar magnetic field and ajcM is the initial pitch 
angle (see [21 ] for more details). Magnetic focusing and magnetic mirroring are the two competing 
effects arising in presence of a magnetic field. The former increases the CR flux where the field 
is more concentrated, enhancing the latter bounces CRs out of the cloud as long as a —> n/l, 
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namely when the CR velocity is perpendicular to the magnetic field line. A simple fitting formula 
which combines in a single expression the effects of energy losses and magnetic fields for different 


CR proton and electron interstellar spectra is given by [124[] : 




„(low N)^(high N) 


.(high N) 
^0,k 


Nia) 


1020 


cm 


^(low N) 
'^^0,k 


exp 


Z(a) 


-0,t: 


-1 


(2.4) 


where 2(a) = jumpA(a) / cos a is the effective surface density seen by a CR propagating with pitch 
angle a, nip the proton mass and p = 2.36 the molecular weight for the assumed fractional abun¬ 
dances of H 2 and He. The fitting coefficients of Eq. (^^ are given in |^, 


3. Cosmic-ray ionisation rate at high densities 

At higher column densities {N > 10^^ cm“2) the drop in becomes even more dramatic, since 
the CR attenuation starts to be exponential. Besides, the stronger the toroidal field component, the 
larger will be the path travelled by a particle. A useful fitting formula to evaluate the effective 
column density covered by a charged particle is given by [|^. If A(H 2 ) is the average column 
density seen by an isotropic flux of CRs, Aeff has the form 

Aeff = (l+27rrOA^(H2). (3.1) 


The factor T depends on the ratio between the toroidal and the poloidal components of the magnetic 
field, h = \B^IBp\, as well as on its module. It reads 


iBi 

IpG 2 


(3.2) 


where 

_ /t 

b* = -(3.3) 

^max ~ ^min 

6min and 6max being the minimum and the maximum value of b in the whole data cube, respectively. 
When the magnetic field strength is negligible, CRs propagate along straight lines. In this case 
T - 0 and - A^(H 2 ), otherwise (F > 0 and Agff > A^(H 2 ). Besides, the higher the density, the 
stronger is the role of the magnetic field in increasing This justifies the presence in Eq. ( pT| ) 
of the power s that reads 

,^0.7 l»gl0(»/»min) ^34^ 

logio(« max /^min) 

?imin and Umax being the minimum and the maximum value of the density in the whole data cube, 
respectively. Notice also that the factor F depends both on the local value of the magnetic field 
(through b and |B|) and on the large scale configuration (by means of b* and s). Once evaluated the 
effective column density, the corresponding (effective) CR ionisation rate, is obtained by 




(3.5) 
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where ^^^(A^eff) is computed using Eq. after replacing N{a) and E(a) with A^eff and yumnA^eff) 
respectively. The factor k is given by 


1 1 /SOOyuG 

K-- + — arctan ——— 

2 ;r \ |B| 


(3.6) 


and it represents the correction for magnetic effects. 


4. Diffusion coefficients and time scales 


Ions and electrons are frozen into magnetic field, while neutrals can diffuse through reaching 
the central part of the core. A frictional force couples charged particles and neutrals. Since CRs 
regulate the ionisation degree, they set limits on the coupling between gas and magnetic field, 
controlling the time scale of the collapse. In general, at lower densities < 10^ - 10^ cm“^) 
the diffusion is controlled by the ambipolar resistivity. At intermediate densities (10^ - 10^ cm 


< 


nii 2 ~ cm“^) Hall diffusion dominates: the more massive charged species (molecular ions and 

charged grains) are coupled with the neutral gas through collisions and then they are decoupled 
from magnetic field. At the highest densities {ny{^ > lO'^ cm“^, 1^]) ions and electrons are knocked 


off from field lines due to collisions and the Ohmic dissipation sets in. However, the extent to 
which a diffusion process dominates over the others hinges on several factors, one of which is the 


assumed grain size distribution (see e.g. Fig. 4 in pSQ) as well as the CR ionisation rate. For this 
reason, after determining the variation of as a function of the density distribution and magnetic 
field configuration, we investigated about how microscopic resistivities are affected by remarkable 
deviation of the CR ionisation rate from a constant value. 

The induction equation reads 


dB ^ ^ ^ B 

— + V X (S X 17) ^ V X 1770V X S + ?7h(V X S) X - + ?7ad 
ot B 


. s 

(VxS)x- 


(4.1) 


where U is the fluid velocity and B the magnetic field vector. Ambipolar, Hall, and Ohmic resis¬ 
tivities (77 ad> ??h> and rjo, respectively) can be written as a function of the parallel ((T||), Pedersen 
(cTp) and Hall (cth) conductivities (e.g. p5[]). These conductivities depend on the ionisation 

fraction that, in turn, is determined by We calculated the ionisation fraction of all the species 
involved (electrons, metal and molecular ions, neutral and charged grains, H^ and Hj) using the 
chemical model described in Appendix A in [|2^. 

The drift velocity of the magnetic field {Ob) can be represented by the velocity of the charged 
species, which are frozen with field lines, with respect to neutrals. From the comparison of this 
velocity with the fluid velocity, it is possible to assess the degree of diffusion of the field and then 
to estimate the size of the region where gas and magnetic field are decoupled. Ub can be written 
as a function of resistivities 0, allowing to isolate the ambipolar (AD), Hall (H), and Ohmic (O) 
contributions, namely 

Ub = 0ad + 0u + Uo, (4.2) 


where 


Uad - 


47r ?7 ad 
cR2 


jxB, 


0n = 


47r?7H 

cB^ 


{jxB)xB, 


Uo = 


4/7 770 V 
cB2 ^ 


xB. 


(4.3) 
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Then 


^ Ttt 


(?7ad + ?7o)/xS + ?7h r 


B 


(4.4) 


where j = {clAn)\ xB is the current density. Thus, the diffusion time of the magnetic field, tg, can 
be written as a function of the time scales associated to the three diffusion processes. 


1 

tB 


1 1 1 

^AD to 


(4.5) 


where tk - RjUk {k - AD,H,0) and R a typical length scale of the region. The diffusion time of 
the magnetic field can then be compared to the time scale of evolution of the fluid, tdyn- We define 
the dynamical time scale of the cloud as t^yn = RIU, where U is the fluid velocity, including both 
infall and rotation. In regions where ts < tdyn the magnetic field is partially decoupled and therefore 
has less influence on the gas dynamics while, if tg > tdynj diffusion is not efficient enough and the 
magnetic field remains well coupled to the gas. 

We performed our calculations for three snapshots from numerical simulations: two low- and 
one high-mass cases. We adopted both a spatially uniform CR ionisation rate of = 5 x 10~'^ s“^, 
and the formalism described in the previous sections to evaluate the attenuation of CRs in a mag¬ 
netised cloud. We found that the ionisation fraction is significantly lower in the second case, and 
therefore the coupling with the magnetic field is weaker than usually assumed in the central region 
of a collapsing cloud (see Fig. 5-7 in [^). 


5. Conclusions 

We applied our modelling for cosmic-ray propagation from large (0.1 pc) to small (< 100 AU) 
scales where the formation of a protostellar disc is expected. We demonstrated that CRs play a 
lead role in regulating the decoupling between the gas and the magnetic field and then in the cloud 
collapse. As an instance, a decoupling zone of radius of 50- 100 AU around the central protostar 
is found in the case with variable but not when is assumed constant for a low-mass case 


snapshot. Besides, this size compares well with the size of the protostellar discs (see Fig. 6 in []25p). 
Also in the high-mass case we examined (see Fig. 7 in [25]) a decoupling zone of size ~ 100 AU is 


found. However, its size becomes smaller for smaller grains, and disappears altogether for grains 
smaller than 10“^ cm. 

Although the models adopted do not represent a time sequence, they nevertheless suggest that 
a decrease in the ionisation and/or an increase in the resistivity occurs in the innermost region of a 
cloud some time after the onset of collapse, but not earlier. In fact, the conditions for a substantial 
increase in the magnetic diffusion time are that the field is considerably twisted and that dust grains 
had time to grow by coagulation. It is tempting to speculate that large, 100 AU-size discs are only 
allowed to form at a later stage when the powerful magnetic brake on the infalling gas has been 
relieved by either (or a combination) of these effects. 

Finally, in addition to the size of the decoupling zone is also determined by the grain 
size assumed in the chemical model and the volume of this region decreases for smaller grain 
size. The maximum grain size assumed in our work (flmin = 10“^ cm) is likely a realistic value 
for the condition expected in disc-forming regions, where larger grains are predicted to form by 
compression and coagulation of smaller grains. 
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